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Passive  direct  methanol  fuel  cells  (DMFCs)  are  promising  energy  sources  for  portable  electronic  devices. 
Different  from  DMFCs  with  active  fuel  feeding  systems,  passive  DMFCs  with  nearly  stagnant  fuel  and  air 
tend  to  bear  comparatively  less  power  densities.  A  steady  state,  one-dimensional,  multi-component  and 
thermal  model  is  described  and  applied  to  simulate  the  operation  of  a  passive  direct  methanol  fuel  cell. 
The  model  takes  into  consideration  the  thermal  and  mass  transfer  effects,  along  with  the  electrochemical 
reactions  occurring  in  the  passive  DMFC.  The  model  can  be  used  to  predict  the  methanol,  oxygen  and 
water  concentration  profiles  in  the  anode,  cathode  and  membrane  as  well  as  to  estimate  the  methanol 
and  water  crossover  and  the  temperature  profile  across  the  cell.  Polarization  curves  are  numerically 
simulated  and  successfully  compared  with  experiments  for  different  methanol  feed  concentrations.  The 
model  predicts  with  accuracy  the  influence  of  the  methanol  feed  concentration  on  the  cell  performance 
and  the  correct  trends  of  the  current  density  and  methanol  feed  concentration,  on  methanol  and  water 
crossover.  The  model  is  rapidly  implemented  and  is  therefore  suitable  for  inclusion  in  real-time  system 
level  DMFC  calculations.  Due  to  its  simplicity  the  model  can  be  used  to  help  seek  for  possibilities  of 
optimizing  the  cell  performance  of  a  passive  DMFC  by  studying  impacts  from  variations  of  the  design 
parameters  such  as  membrane  thickness,  catalyst  loading,  diffusion  layers  type  and  thicknesses. 
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1.  Introduction 

Conventional  batteries  are  becoming  inadequate  for  the  increas¬ 
ing  power  requirements  of  portable  electronic  devices  such  as 
mobile  phones,  PDA’s,  laptops  and  multimedia  equipment.  Direct 
methanol  fuel  cells  (DMFCs)  are  promising  candidates  as  portable 
power  sources  because  they  do  not  require  any  fuel  processing 
resulting  in  a  simpler  design  and  operation,  higher  reliabil¬ 
ity  and  operate  at  low  temperatures.  DMFCs  offer  high  energy 
densities,  longer  runtime,  instant  recharging  and  lower  weight 
than  conventional  batteries.  The  most  significant  obstacle  for 
DMFC  development  is  methanol  crossover,  since  methanol  diffuses 
through  the  membrane  generating  heat  but  no  power.  This  problem 
can  be  limited  if  the  cell  operates  with  low  methanol  concentra¬ 
tion  on  the  anode.  Flowever,  this  significantly  reduces  the  system 
energy  density  since  water  will  produce  no  power  and  will  take 
up  a  large  volume  in  the  fuel  reservoir.  Due  to  the  concentration 
gradient  between  the  anode  and  cathode,  water  crosses  through 
the  membrane.  The  presence  of  a  large  amount  of  water  floods  the 
cathode  and  reduces  cell  performance. 
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There  are  two  types  of  fuel  and  oxidant  supply  in  a  DMFC:  an 
active  and  a  passive  one.  Active  systems  use  extra  components  such 
as  a  pump  or  blower,  a  fan  for  cooling,  reactant  and  product  con¬ 
trol,  which  allows  the  operation  of  a  DMFC  at  favourable  conditions 
with  respect  to  temperature,  pressure,  concentration  and  flow  rate. 
This  type  of  system  supply  is  more  complex,  has  greater  costs  and 
lower  system  energy  densities.  More  recently,  the  passive  DMFCs 
have  been  proposed  and  investigated  [1-17].  Passive  systems  use 
natural  capillary  forces,  diffusion,  convection  (air  breathing)  and 
evaporation  to  achieve  all  processes  without  any  additional  power 
consumption.  Therefore,  the  fuel  cell  system  becomes  much  sim¬ 
pler  and  more  compact  being  more  suitable  for  portable  power 
sources.  The  passive  DMFCs  have  much  lower  power  density  due 
to  the  inability  to  handle  the  excess  water  produced  on  the  cathode 
and  crossed  from  the  anode  and  to  the  excess  of  heat  lost  from  the 
fuel  cell  to  the  ambient  air.  Therefore,  the  key  issues  in  the  portable 
DMFC  system  is  the  thermal  and  water  management  [5,12,14-17]. 

Performance  of  a  passive  DMFC  relies  on  a  vast  number  of 
parameters,  including  the  methanol  feed  concentration,  efficien¬ 
cies  of  methanol  and  oxygen  transport  within  the  different  layers, 
the  release  rate  of  gaseous  carbon  dioxide  and  its  effect  on  methanol 
transport,  the  specific  area  of  catalyst  in  the  catalyst  layers,  the 
thickness  of  the  membrane,  the  gas  diffusion  layer  properties,  the 
rate  of  methanol  and  water  permeation  and  so  on.  Experimental 
investigating  of  the  impact  of  these  parameters  one  by  one  through 
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Nomenclature 

a 

specific  surface  area  of  the  anode,  cm-1 

Aq 

active  area,  cm2 

A, 

total  area  without  the  holes,  cm2 

Wholes 

total  area  of  the  holes,  cm2 

c 

concentration,  mol  cm-3 

c2 

concentration  at  the  AAP/ACP  interface,  mol  cm-3 

C3 

concentration  at  the  ACP/AD  interface,  mol  cm-3 

C4 

concentration  at  the  AD /AC  interface,  mol  cm-3 

C5 

concentration  at  the  AC/membrane  interface, 
mol  cm-3 

C6 

concentration  at  the  membran e/CC  interface, 
mol  cm-3 

C7 

concentration  at  the  CC/CD  interface,  mol  cm-3 

C8 

concentration  at  the  CD/CCP  interface,  mol  cm-3 

c0  2,  ref 

reference  concentration  of  oxygen,  mol  cm-3 

cP 

specific  heat  capacity,  J  mol-1 1<-1 

dE/dT 

rate  of  change  of  electromotive  force,  VK-1 

D 

diffusion  coefficient,  cm2  s-1 

DeJ? 

effective  diffusion  coefficient,  cm2  s-1 

Ecell 

thermodynamic  equilibrium  potential,  V 

F 

FaradayTs  constant,  96,500  C  mol-1 

G 

Gibbs  free  energy,  J  mol-1 

g 

gravitational  acceleration,  cm2  s-1 

H 

enthalpy  of  reaction,  J  mol-1 

hmass 

mass  transfer  coefficient,  cms-1 

hheat 

heat  transfer  coefficient,  W  cm-2  K-1 

kell 

cell  current  density,  A  cm-2 

Ich3oh 

leakage  current  density  due  to  methanol  crossover, 

A  cm-2 

tCH3  OH 
0,ref 

exchange  current  density  of  methanol,  A  cm-2 

/°2 

0,ref 

exchange  current  density  of  oxygen,  A  cm-2 

J'a 

volumetric  current  density,  A  cm-3 

k 

constant  in  the  rate  expression  (Eq.  (27)) 

& 

1 

00 

partition  coefficients 

I< 

thermal  conductivity,  W  cm-1  K-1 

L 

length  of  the  active  area,  cm 

nd 

electro-osmotic  drag  coefficient  of  water 

N 

molar  flux,  mol  cm-2  s-1 

Pair 

pressure  of  air  in  cathode,  atm 

qAC 

heat  generated  in  AC,  W  cm-2 

Qcc 

heat  generated  in  CC,  W  cm-2 

Q 

heat  transfer,  W 

R 

gas  constant,  8.314Jmol-1  K-1 

P-Cell 

internal  resistance  of  the  fuel  cell,  cm2  S-1 

Pcond 

conduction  resistance,  K  W_1 

Pconv 

convection  resistance,  KW-1 

Ptotal 

total  thermal  resistance,  KW-1 

T 

temperature,  K 

Uch3oh 

thermodynamic  equilibrium  potential  of  methanol 
oxidation,  V 

U02 

thermodynamic  equilibrium  potential  of  oxygen 
oxidation,  V 

Vcell 

cell  voltage,  V 

*ch3oh 

mole  fraction  of  methanol,  mol  mol-1 

X 

coordinate  direction  normal  to  the  anode,  cm 

Greek 

A 

variation 

a 

net  water  transport  coefficient 

aA 

anodic  transfer  coefficient 

ac 

cathodic  transfer  coefficient 

P 

coefficient  of  volume  expansion,  1  K-1 

8 

thickness,  cm 

£ 

porosity 

1 

overpotential,  V 

K 

ionic  conductivity  of  the  membrane,  S  cm-1 

X 

constant  in  the  rate  expression  (Eq.  (27)),  mol  cm-3 

fi 

dynamic  viscosity,  gem-1  s-1 

V 

kinematic  viscosity,  cm2  s-1 

P 

density,  gem-3 

Vo2 

stoichiometric  coefficients  of  oxygen  in  the  cathode 
reaction 

vH20 

stoichiometric  coefficients  of  water  in  the  cathode 
reaction 

vcross,02 

stoichiometric  coefficients  of  oxygen  in  the  unde¬ 
sired  cathode  reaction 

V  cross, H2C 

,  stoichiometric  coefficients  of  water  in  the  unde¬ 
sired  cathode  reaction 

$ch3oh 

electro-osmotic  drag  coefficient  of  methanol 

Subscripts 

A 

anode 

air 

air 

C 

cathode 

CH3OH 

methanol 

i 

species  i 

j 

species  j 

H20 

water 

02 

oxygen 

Superscripts 

0 

feed  conditions 

AAP 

anode  acrylic  plate 

ACP 

anode  copper  plate 

AC 

anode  catalyst  layer 

AD 

anode  diffusion  layer 

CAP 

cathode  acrylic  plate 

CC 

cathode  catalyst  layer 

CCP 

cathode  copper  plate 

CD 

cathode  diffusion  plate 

l 

plate  l 

M 

membrane 

t 

plate  t 

is  not  time  or  cost  efficient.  In  order  to  help  understand  the  oper¬ 
ation  of  a  passive  DMFC  and  locate  the  key  parameters  on  cell 
performance,  a  theoretical  model  is  essential. 

In  a  previous  work  Oliveira  et  al.  [18]  report  an  intensive  review 
on  the  work  done  in  DMFC  empirical  and  fundamental  modelling. 
Despite  the  number  of  studies  in  DMFCs  modelling  only  a  few  sim¬ 
ulate  passive  DMFCs  [6,19-23]  and  only  a  small  part  took  into 
account  thermal  effects  [6,19,22]. 

Since  thermal  management  is  a  key  issue  in  the  portable  DMFC 
system  it  is  important  to  develop  new  models  accounting  for  this 
effect  and  that  can  be  a  simple  computer-aided  tool  to  the  design 
and  optimization  of  passive  direct  methanol  fuel  cells. 

Chen  and  Zhao  [19]  presented  a  one-dimensional  model  to 
describe  a  passive  liquid-feed  direct  methanol  fuel  cell  combin¬ 
ing  the  effects  of  heat  and  mass  transfer.  The  model  provides  the 
temperature  profile  along  the  different  layers  of  a  passive  DMFC. 

More  recently,  Chen  et  al.  [22]  presented  a  two-dimensional 
two-phase  thermal  model  for  passive  direct  methanol  fuel  cells 
(DMFC).  The  model  was  based  on  the  unsaturated  flow  the¬ 
ory  in  porous  media.  The  model  is  solved  numerically  using  a 
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Fig.  1.  Schematic  representation  of  a  passive  DMFC. 


home-written  computer  code  to  investigate  the  effects  of  various 
operating  and  geometric  design  parameters,  including  methanol 
concentration  as  well  as  the  open  ratio  and  channel  and  rib  width 
of  the  current  collectors,  on  cell  performance. 

As  disadvantages,  these  two  models  [19,22]  consider  the  cat¬ 
alyst  layers  as  an  interface,  so  it  is  not  possible  to  obtain  the 
temperature  and  concentration  profiles  in  these  layers,  and  the 
authors  assumed  that  the  anode  side  is  well  insulated  so  no  heat  is 
lost  from  the  anode  side.  This  assumption  may  be  very  unrealistic 
in  a  passive  DMFC  working  in  a  portable  system. 

Based  on  the  model  developed  previously  by  the  same  authors 
[24],  the  goal  of  the  present  work  is  the  development  of  a  steady 
state,  one-dimensional,  multi-component  and  thermal  model.  The 
model  takes  into  consideration  the  thermal  and  mass  transfer 
effects,  along  with  the  electrochemical  reactions  occurring  in  the 
passive  DMFC.  The  model  can  be  used  to  predict  the  methanol,  oxy¬ 
gen,  carbon  dioxide  and  water  concentration  profiles  in  the  anode, 
cathode  and  membrane  as  well  as  to  estimate  the  methanol  and 
water  crossover  and  the  temperature  profile  across  the  cell.  The  aim 
of  the  work  is  to  produce  a  simplified  model  describing  the  main 
heat  and  mass  transfer  effects  in  a  passive  DMFC  fuel  cell  and  repro¬ 
ducing  with  satisfactory  accuracy  experimental  data.  The  results  of 
a  simulation  study  using  a  developed  model  for  passive  DMFCTs  are 
presented.  The  model  was  validated  with  data  from  experiments 
conducted  in  an  in-house  designed  passive  DMFC  and  with  recent 
published  data  [3]. 

2.  Model  development 

A  schematic  representation  of  a  passive-feed  direct  methanol 
fuel  cell  is  shown  in  Fig.  1,  consisting  of 

•  an  acrylic  plate  (AAP)  containing  the  fuel  tank,  a  copper  plate 
(ACP),  a  diffusion  layer  (AD)  and  a  catalyst  layer  (AC)  at  the  anode 
side; 

•  a  polymer  electrolyte  membrane  (M); 

•  a  catalyst  layer  (CC),  a  diffusion  layer  (CD),  a  copper  plate  (CCP), 
and  an  acrylic  plate  (CAP)  at  the  cathode  side. 


The  direct  methanol  fuel  cell  is  complex  system  involving  simul¬ 
taneous  mass,  charge  and  energy  transfer.  In  order  to  simplify  the 

processes  occurring  in  a  DMFC  the  following  simplifications  and 

assumptions  were  made: 

•  the  fuel  cell  is  assumed  to  operate  under  steady-state  conditions; 

•  the  transport  of  heat  and  mass  through  the  gas  diffusion  and  cat¬ 
alyst  layers  is  assumed  to  be  a  diffusion-predominated  process 
and  the  convection  effect  is  negligible; 

•  mass  transport  in  the  diffusion  layers  and  membrane  is  described 
using  effective  Fick  models; 

•  the  thermal  energy  model  is  based  on  the  differential  thermal 
energy  conservation  equation  (Fourier’s  law); 

•  pressure  gradient  across  the  layers  is  negligible; 

•  only  the  liquid  phase  is  considered  in  the  anode  side,  so  carbon 
dioxide  remains  dissolved  in  solution; 

•  gaseous  methanol  and  water  are  considered  in  the  cathode; 

•  solutions  are  considered  ideal  and  dilute; 

•  local  equilibrium  at  interfaces  is  represented  by  partition  func¬ 
tions; 

•  the  catalyst  layers  are  assumed  to  be  a  macro-homogeneous 
porous  electrode  so  reactions  in  these  layers  are  modelled  as  a 
homogeneous  reaction; 

•  anode  kinetics  is  described  by  step  mechanism,  with  a  rate 
expression  similar  to  the  used  by  Meyers  and  Newman  [25]; 

•  the  anodic  and  cathodic  overpotential  is  constant  through  the 
catalyst  layers; 

•  cathode  kinetics  is  described  by  Tafel  equation; 

•  methanol  and  water  transport  through  the  membrane  is  assumed 
to  be  due  to  the  combined  effect  of  the  concentration  gradi¬ 
ent  between  the  anode  and  the  cathode  and  the  electro-osmosis 
force; 

•  on  the  anode  side,  the  heat  and  mass  transfer  of  methanol  from 
the  bulk  solution  to  the  ACP  is  assumed  to  be  driven  by  natural 
convection; 

•  on  the  cathode  side,  the  heat  and  mass  transfer  between  the  CCP 
and  the  ambient  occur  by  natural  convection; 

•  the  heat  generation  by  electrochemical  reactions  occurring  in  the 
catalyst  layers  is  considered; 

•  when  compared  with  the  heat  generated  by  electrochemical  reac¬ 
tions  and  overpotential,  the  heat  released  by  joule  effects  is 
ignored; 

•  the  temperatures  of  the  external  walls  of  the  cell  (T0  and  T9  in 
Fig.  1)  are  known; 

•  the  heat  flux  generated  in  the  catalyst  layers  is  assumed  to  be 
constant. 

2.1.  Mass  transport 
Anode  reaction: 

Methanol  oxidation  :  CFI3OFI  +  FI20  C02  +  6H+  +  6e_ 


In  a  passive-feed  DMFC  the  fuel,  methanol  or  an  aqueous 
methanol  solution,  and  the  oxidant  are  supplied  to  the  reaction 
zone  by  natural  convection.  From  the  ACP  through  the  AD  and  from 
the  AC  through  the  M,  methanol  solution  is  transported  primarily 
by  diffusion.  In  a  similar  way  the  transport  of  oxygen  on  the  CCP, 
CD  and  CC  is  enhanced  by  diffusion.  After  the  electrochemical  reac¬ 
tion  of  methanol  oxidation,  witch  takes  place  in  the  AC,  the  carbon 
dioxide  produced  moves  counter-currently  toward  the  AAP.  At  suf¬ 
ficiently  high  current  densities  carbon  dioxide  emerges  in  the  form 
of  gas  bubbles  from  the  surface  of  the  AC.  In  the  CC,  oxygen  reacts 
with  protons  and  electrons  generating  water.  The  water  produced 
in  CC  moves  counter-currently  toward  the  CCP  and  also  under  some 
operating  conditions,  by  back  diffusion  toward  the  anode. 


Cathode  reaction: 

Oxygen  reduction  :  02  +  4H+  +  4e_  ->  2FI20 


Methanol  oxidation  :  CFI3OFI  +  (3/2)02  ->  C02  +  2H20 


The  transport  process  of  methanol  and  water  from  the  fuel  tank 
to  the  ACP  are  described  by 


N-  _  JjAAP  _  CAAP\ 

~  nmass,j^j  S  > 


(1) 


where  j  represents  methanol  or  water,  N  is  the  molar  flux,  C  is  the 
molar  concentration  and  hmass  is  the  mass  transfer  coefficient. 

In  the  anode  copper  plate,  diffusion  and  catalyst  layer,  the 
methanol  and  water  flux  are  related  to  the  concentration  gradi- 
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ent  by  assuming  Fickian  diffusion  [26]  with  an  effective  diffusivity 
DefACP  in  the  ACP,  D'ffAD  in  the  AD  and  Df'AC  in  the  AC.  The 
methanol  and  water  flux  can  be  determined  from: 

dCACP 

Nj  =  -Dj  — ,  j  represents  methanol  or  water,  (2) 

dCfD 

Nj  =  -DjfL  j  represents  methanol  or  water  (3) 

and 

dCAC 

Nj  =  -Djff,AC  »  J  represents  methanol  or  water  (4) 

The  concentration  at  the  AAP/ACP,  ACP/ AD  and  AD /AC  interfaces 
is  given  by  assuming  local  equilibrium  with  a  partition  coefficient 
I<2,I<3  and  I<4 ,  respectively.  The  boundary  conditions  for  Eqs.  (2)-(4) 
are  (see  Fig.  1) 


Atx  =  x2  :  CA cp  =  I<2Caap,  j  represents  methanol  or  water  (5) 

Atx  =  x3  :  CADj  =  K3Cacp  ,  j  represents  methanol  or  water  (6) 

Atx  =  x4  :  Cfj  =  I<4CADj,  j  represents  methanol  or  water  (7) 

Atx  =  X5  :  CAC  =  CACj ,  j  represents  methanol  or  water  (8) 

In  fuel  cells,  all  the  fluxes  can  be  related  to  a  single  character¬ 
istic  flux,  the  current  density  or  charge  flux  of  the  fuel  cell.  In  the 
DMFC,  the  methanol  flux  is  related  to  the  current  density  and  the 
permeation  flux  of  methanol  through  the  membrane,  (N^30H),  by: 

nch3oh  =  (9) 

where  F  represents  the  Faraday’s  constant  and  /ceH  is  the  cell  current 
density. 

At  the  anode  side,  the  water  flux  is  related  to  the  current  density 
and  to  the  net  water  transport  coefficient,  a  (defined  as  the  ratio 
of  the  net  water  flux  though  the  membrane  from  the  anode  to  the 
cathode  normalized  by  protonic  flux),  by: 

%2o  =  !f(a  +  l)  (10) 


The  transport  of  methanol  and  water  through  the  membrane  is 
assumed  to  be  due  to  the  combined  effect  of  the  concentration  gra¬ 
dient  and  the  electro-osmosis  force.  The  fluxes  can  be  determined 
from: 


NM  _  _Deff,M 

nu  —  ru  n 


dC™ 


yCH3OH 


CH3OH  t  ICell 
+?Of3GH-r 


1  dCM  1 

NUn  =a!m  =  +  nd  !m 


6  F 


H2 0  dX 


(11) 

(12) 


diffusion  [26]  with  an  effective  diffusivity  Dy^,cc.  The  flux  can  be 
determined  from: 

dCcc 

Nj  =  -Djff,cc  j  represents  methanol,  water  or  oxygen  (14) 

It  is  here  considered  that  the  entire  methanol  crossing  the  mem¬ 
brane  reacts  at  the  cathode  catalyst  layer  so  the  concentration  at 
the  CC/CD  interface  is  zero.  It  is  assumed  that  there  is  no  oxygen 
crossover,  so  the  oxygen  concentration  in  CC/M  interface  is  zero. 
The  concentration  of  water  and  methanol  at  the  membrane/CC 
interface  and  the  concentration  of  water  and  oxygen  at  the  CC/CD 
interface  are  given  by  assuming  local  equilibrium  with  a  partition 
coefficient  I<6  and  JC7,  respectively.  The  boundary  conditions  for  Eq. 
(14)  are: 

Atx  =  x6  :  Cfj  =  I<6 C^j  j  represents  methanol  or  water  and 

xC6CCo2=  0  (15) 

Atx  =  x-j  :  Cga0H  =  0,  C^0  =  CjCH2q  and  (16) 

At  the  cathode  catalyst  layer,  the  oxygen  reacts  with  the  elec¬ 
trons  and  protons  to  produce  water.  However,  part  of  oxygen  fed 
is  consumed  due  to  methanol  crossover  to  form  an  internal  current 
and  a  mixed  potential.  Therefore  the  oxygen  flux  is  related  to  the 
current  density  and  the  permeation  flux  of  methanol  through  the 
membrane  by: 

N02  =  Vo2  +  vcross,02^CH3OH  (1^) 

where 

V02  —  1  cUld  O  cross,  02  — 

The  vq2  represents  the  stoichiometric  coefficient  of  oxygen  in  the 
cathode  reaction  and  the  vcrosSjo2  stoichiometric  coefficient  of  oxy¬ 
gen  in  the  undesired  cathode  reaction. 

At  the  cathode  side,  the  water  flux  is  related  to  the  water  produc¬ 
tion  from  the  oxygen  reduction  reaction  and  methanol  crossover 
oxidation  and  to  the  net  water  flux  transported  from  the  anode  to 
the  cathode  by: 

Nh20  =  vH2o^  +  vcr oss,H2oNchs0h  +  K2o  (18) 

where 

vcross,H20  =  2,  Vh2q  =  2 

The  vH2o  represents  the  stoichiometric  coefficient  of  water  in  the 
cathode  reaction  and  the  vcwss  H2o  stoichiometric  coefficient  of 
water  in  the  undesired  cathode  reaction. 

In  the  cathode  diffusion  layer  and  cathode  copper  plate  the  oxy¬ 
gen  and  water  flux  are  related  to  the  concentration  gradient  by 


The  electro-osmotic  drag  ($ch3oh,  in  Eqs.  (11)  and  (12),  is 
defined  as  the  number  of  methanol  or  water  molecules  dragged  by 
the  hydrogen  ions  moving  through  the  membrane. 

The  net  water  transport  coefficient,  a,  can  be  calculated  using 

Eq*  (12). 

The  concentration  at  the  AC/membrane  interface  is  given  by 
assuming  local  equilibrium  with  a  partition  coefficient  I< 5.  The 
boundary  conditions  for  the  integration  of  equations  Eqs.  (11)  and 
(12)  is  given  by 

Atx  =  x5  :  fj  =  K3CACj  j  represents  methanol  or  water  (13) 

In  the  cathode  catalyst  layer,  the  methanol,  water  and  oxygen 
flux  are  related  to  the  concentration  gradient  by  assuming  Fickian 


JV,.  =  -DfCD 


dC[D 

dx 


N,  =  -DfCCP 


dCfcp 

dx 


represents  oxygen  or  water  vapour 
i  represents  oxygen  or  water  vapour 


(19) 

(20) 


where  Df^,CD  and  Df^,ccp  are  the  effective  diffusion  coefficient  of 
oxygen  and  water  in  the  CD  and  CCP. 

The  concentration  at  the  CC/CD  and  CD/CCP  interfaces  is  given 
by  assuming  local  equilibrium  with  a  partition  coefficient  I<7  and 
I<8.  The  boundary  conditions  for  Eqs.  (19)  and  (20)  are: 


At  x  =  x7  :  C™  =  Kj  jCjC  i  represents  oxygen  or  water  vapour 

(21) 
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At  x  =  Xg  :  C8c7  =  Kg  jCgD  i  represents  oxygen  or  water  vapour 

(22) 


Like  at  the  anode  side,  the  transport  process  of  oxygen  from  the 
air  to  the  CCP  is  described  by: 


N- —  hc  (C°-Cccp) 
iV*  “  “mass, AS'  l9 ,i  ' 


(23) 


where  i  represents  oxygen. 

We  assume  that  the  air  at  the  CCP  is  in  a  saturated  state,  then  the 
water  vapour  feed  concentration  (C°  H20)  is  equal  to  water  vapour 

concentration  (Cg,s^0)  and  can  be  determined  from  the  saturated 
pressure  or  moist  air. 

To  account  for  the  effect  of  methanol  crossover  on  the  cathode 
overpotential  it  is  assumed  that  the  methanol  crossing  the  mem¬ 
brane  completely  reacts  electrochemically  at  the  cathode.  In  this 
way  the  internal  current  (/ch3oh)  due  to  methanol  oxidation  can  be 
written  as: 


Ich3oh  =  S™ch3oh  (24) 

where  the  methanol  flux  in  the  membrane  (N^H  0H)  is  obtained 
from  Eq.  (11). 

The  volumetric  current  density  (jA )  expression  for  methanol  oxi¬ 
dation  is  taken  from  Meyers  and  Newman  [25]  as 


Ja 


kCAC 

aICH3OH _ /CLCH3OH _ 

°’ref  cCH3OH  +  *  expia atjaF I RTac) 


(25) 


where  a  represents  the  specific  surface  area  of  the  anode,  1^/ H  *s 
the  exchange  current  density  of  methanol,  aA  is  the  anodic  transfer 
coefficient,  rjA  is  the  anode  overpotential,  TAC  is  anode  catalyst  layer 
temperature,  and  k  and  X  are  constants. 

The  current  density  is  related  to  the  volumetric  current  density 
using  the  following  equation 


kell  = 


r*  e  r*  e 

=  Ja  =  /  < 

Jx5  Jx5 


afH3OH _ 

ai0  ,ref  c/KC 


kCAC 

CH3OH 


+  X  exp (aArjAF/RTAC) 


ch3oh 


x  exp 


(26) 


Eq.  (26)  is  used  to  calculate  the  anode  overpotential  for  a  given  Icdh 
assuming  rjA  as  constant  in  the  anode  catalyst  layer  AC. 

At  the  cathode,  the  electrochemical  reaction  is  modelled  using 
Tafel  equation  for  the  oxygen  reduction  taking  in  account  the  mixed 
potential.  The  cathode  overpotential  can  then  be  determined  from: 


Icell+lcH3OH  =  I^ref^~  exP  (^~)  (27) 

02,  ref  K  LL 

where  1° 2.^  represents  the  exchange  current  density  of  oxygen,  ac 
is  the  cathodic  transfer  coefficient,  pc  is  the  cathode  overpotential 
and  TCc  is  the  cathode  catalyst  layer  temperature. 

The  mass  transfer  coefficient  in  Eqs.  (1)  and  (23)  can  be  deter¬ 
mined  from  [27]: 


Sh  = 


D 


0.825  + 


0.387  x  Ra h6 


(1  +  (0.492/Sc  )9/16)8/27 


(28) 


where  Ra  is  the  Rayleigh  number  (Ra  =  Gr  x  Sc),  Sc  is  the  Schmidt 
number  (Sc  =  v/D),  Gr  is  the  Grashof  number  (Gr =gA CL3 /Cu2 ),  L  rep¬ 
resents  the  length  of  the  active  area,  D  is  the  diffusion  coefficient, 
g  is  the  gravitational  acceleration  and  v  is  the  kinematic  viscosity. 


2.2.  Heat  transport 


Based  on  the  simplifications  and  assumptions  described  previ¬ 
ously  the  following  overall  heat  transfer  equation  can  be  proposed 
(see  Fig.  1): 

0.AC  +  Qcc  =  Q,  +  02  (29) 

The  total  heat  generated  in  the  DMFC  is  equal  to  the  heat  losses 
to  the  surrounding  environment  at  the  anode  and  cathode. 

Complementarily,  the  following  heat  transfer  balances  can  be 
written: 

Q3  =  Q.ac  -  Qi  (30) 

Q2=Qcc  +  Q3  (31) 

At  the  anode,  heat  generated  by  the  electrochemical  reaction  in 
the  AC  is  given  by 

Qac  =  Iceim  -  ken  (AHAqFAGa)  (32) 

In  this  equation  the  first  term  represents  the  heat  due  to  the 
activation  and  mass  transfer  overpotentials  at  the  anode  and  the 
second  term  represents  the  entropy  change  of  the  anodic  electro¬ 
chemical  reaction,  with  A HA  denoting  the  anodic  reaction  enthalpy 
and  A GA  the  Gibbs  free  energy. 

In  a  similar  way,  the  heat  generated  at  the  CC,  can  be  determined 
from 

QCC  =  Ucell  +  IcH3Oh)PC  -  (I Cell  +  ^CH3Oh)  ^ ^ 

-Ich3oh  (  AHA6FACa  )  (33) 

where  the  first  term  represents  the  heat  due  to  the  activation 
and  mass  transfer  overpotentials  and  mixed  potential  caused  by 
methanol  crossover  through  the  cathode  and  the  second  term 
represents  the  entropy  change  of  the  cathodic  electrochemical 
reaction,  with  A Hc  denoting  the  cathodic  reaction  enthalpy  and 
A Gc,  the  Gibbs  free  energy  and  the  third  term  denotes  the  entropy 
change  of  methanol  oxidation  reaction  on  the  cathode  due  to 
methanol  crossover. 

In  the  anode  acrylic  plate  section  I  and  diffusion  layer  the  heat 
flux  Qi  can  be  related  to  the  temperature  gradient  across  each  layer, 
using  the  Fourier’s  law,  as 

<34) 

where  l  represents  AAPsectionI  or  AD,  K  is  the  thermal  conductivity 
and  Aa  represents  the  active  area. 

In  the  anode  acrylic  plate  section  II  the  heat  flux  Qi  can,  also,  be 
related  to  the  temperature  gradient  across  this  layer,  using  New¬ 
ton’s  law,  as 


Q  =  -hheatAaAT  (35) 

where  hheat  represents  the  heat  transfer  coefficient. 

At  the  cathode  side  and  membrane,  the  heat  fluxes  Q2  and  Q3 
can  be  related  to  the  temperature  gradient  across  the  CD  and  M 
layers  as 

Q  =  -/(%,  ^  (36) 

where  t  represents  CD  or  M. 

In  a  passive  DMFC  the  copper  plate  has  holes  machined  on  the 
surface,  to  allow  the  reactant  to  reach  the  catalyst  layers  (Fig.  1). 
The  establishment  of  the  heat  transport  equations,  in  this  layer, 
involved  the  consideration  of  two  zones.  In  one  zone  the  heat  is 
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Table  1 

Parameter  values. 


Parameter 

Value 

Reference 

Uo2 

1.24  V 

[28] 

Uch2oh 

0.03  V 

[28] 

dE/dT 

—  1.4  x  10-4  VK-1 

[19] 

K 

0.0365  cm-1 

[28] 

8M 

0.018  cm 

[28] 

foAAPsectionl  fiAAPsectionll 

0.50  cm 

Assumed 

8AD,8CD 

0.015  cm 

Assumed 

SAC,SCC 

0.0023  cm 

Assumed 

sAD,sCD 

0.71 

Assumed 

e* c 

0.81 

Assumed 

ecc 

0.86 

Assumed 

a 

1000  cm-1 

[28] 

tCH3OH 

0,ref 

9.425  x  10-3  exp((35,  570/R)(l/353  -  1/7))  Acm”2 

[19] 

\02 

0  ,ref 

4.222  x  10-6  exp((73,  200/R)(  1/353  -  l/7))Acirr2 

[19] 

k 

7.5  x  10-4 

[28] 

X 

2.8  x  10“9molcm“3 

[28] 

OLa 

0.52 

[28] 

OLc 

1.55 

[28] 

K 2-5 

0.8 

Assumed 

1^7-8, On 

1.25 

Assumed 

K 6 

0.001 

Assumed 

I<7-8,H20 

0.8 

Assumed 

L 

5  cm 

Assumed 

neff,CD,CC 

°2 

neff,CCP 

°2 

£CD,CC2-5[(T 1.75  x  5  8  x  10_4)/(27.772  x  P)]cm2  s”1 

[30] 

[(71-75  x  5.8  x  10-4)/(27.772  x  P)]  cm2  s~3 

[30] 

D^cp 

CH30H 

[(7.608  x  10“7  x  T)I(iiH2o  x  9.485)]  cm2  s”1 

[30] 

neff,AD,AC 
ch3  oh 

^D^c2-5  [(7.608  x  10~7  x  T)l(fiH2o  x  9.485)]  cm2  s”1 

[30] 

Deff,CC 

CH3OH 

ecc 2  5  [(T1-75  x  5.8  x  10“4)/(33.904  x  P)]  cm2  s”1 

[30] 

neff,M 

CH3OH 

4.9  x  10-6exp(2436  x  (1/333  -  1/7))  cm2  s”1 

[28] 

D^CP 

h2o 

[(6.295  x  10-7  x  T)I(/jlCh3oh  x  5.833)]  cm2  s”1 

[30] 

neff,AD,AC 

h2o 

neff,CD,CC 

h2° 

H20 

sad,ac2- 5  [(6.295  x  10~7  x  7)/(/xCh3oh  x  5.833)]  cm2  s”1 

[30] 

£cd,cc2-5[(t i.75  x  6  2  x  10_4)/(25.523  x  P)]cm2  s”1 

[30] 

2.0  x  10-6  exp(2060  x  (1/303  -  1/7))  cm2  s”1 

[30] 

%ch3oh 

2.5  x  Xch3oh 

[28] 

nd 

2.9  exp(1029  x  (1/333  -  1/7)) 

[19] 

fiACP  fiCCP 

0.05  cm 

Assumed 

I<M  ’ 

0.0043  W  cm-1  K_1 

[29] 

kad 

1.95  +  6.57  x  10-47Wm“1  K1 

[29] 

KCD 

1.71+2.96  x  10-57Wm-1  K_1 

[29] 

I<ac 

(1  -eAC)  x  86.7 +  £^(0.341  +  9.26  x  10“4)Wm-1  K"1 

[29] 

K* 

(1  -ecc)x71+£cc(0.0034  +  7.60x  lO-^Wm-1  K_1 

[29] 

transferred  by  conduction  and  in  the  other  (holes)  the  heat  is  trans¬ 
ferred  by  convection.  Using  the  thermal  resistance  concept  [27]  we 
get: 


Q  = 


AT 

Rtotal 


(37) 


where 


The  differential  equations  describing  the  temperature  profiles 
in  the  anode  and  cathode  catalyst  layers  are: 


d2T  _  (lAC 
dx2  I<AC8AC 

d2T  _  Qcc 
dx 2  ~  KCC8CC 


(42) 


(43) 


1 


^total 
Rcond  — 


1  1 

-  +  7^— 


Rcond 

8 

A1  x/C 


since  the  resistances  are  in  parallel  (38) 


(39) 


where  Q/iC  and  Qcc  are,  respectively,  the  heat  generated  in  the  anode 
catalyst  layer  and  cathode  catalyst  layer. 

The  boundary  conditions  for  Eqs.  (42)  and  (43)  are  the  temper¬ 
atures  at  the  walls  (T4,  T5,  T6  and  T7). 

For  these  layers.  Fourier’s  law  gives 


Rconv  — 


Wholes  x  hheat 


and 


(40) 


Atx  =  x4  :  Qi  =  -KACAa  ‘f- 
dx 

rr  dT 

Atx  =  x6  :  Q3  =  -KccAa-r 


(44) 

(45) 


Aa  — A ]  +  A/,0;es  (41) 

where  Rtotai  represents  the  total  thermal  resistance,  Rcond  is  the 
conduction  resistance,  RCOnv  is  the  convection  resistance,  8  is  the 
thickness,  Ahoies  is  the  total  area  of  the  holes  and  Ai  is  the  total  area 
without  the  holes. 


where  dT/dx  is  calculated  using  the  temperature  profile  obtained 
from  the  integration  of  Eqs.  (42)  and  (43). 

Finally,  the  heat  transfer  from  the  AAPsection  I  and  CCP  to  the 
ambient  air  can  be  described  using  the  Newton’s  cooling  law  as 

Q.  =  —hfreatAaAT 


(46) 


V.B.  Oliveira  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  8973-8982 


8979 


The  heat  transfer  coefficient,  due  to  natural  convection  in  Eq. 
(35),  can  be  determined  from  [27]: 


Nu  = 


riheatL 

I< 


0.825  + 


0.387  x  Ra1/6 


(1  +  (0.492/PT)9 


(47) 


where  Ra  is  the  Rayleigh  number  ( Ra  =  Gr  x  Pr),  Pr  is  the  Prandtl 
number  (Pr=  v/K)  and  Gr  is  the  Grashof  number  (Gr=gfiATL3 /v2). 


2.3.  Cell  performance 


The  determination  of  methanol  and  oxygen  concentrations  at 
the  catalyst  layers,  the  temperature  profiles  and  the  anodic  and 
cathodic  overpotentials  from  the  model  equations  enables  predic¬ 
tion  of  the  cell  voltage,  which  can  be  expressed  as: 

VCell  =  ECell  -  Pa  - Pc-  kellRCell  (48) 

where 

Ecell  =  Uo2  -  Uch3oh  +  A T  (Jpj  (49) 

pA  and  pc  are  the  anode  and  cathode  overpotencials  and  the  mem¬ 
brane  resistance  RCeu  is  given  by 

SM 

RCell  =  —  (50) 

where  8M  is  the  membrane  thickness  and  k  is  the  ionic  conductivity 
of  the  membrane. 

The  results  presented  in  the  next  section  were  obtained  based 
on  the  parameters  listed  in  Table  1. 

3.  Experiment 


The  experimental  fuel  cell  consists  of  two  acrylic  end  plates 
(open  on  the  cathode  side  and  with  a  reservoir  on  the  anode  side), 
two  isolating  plates,  two  gold  plated  copper  connector  plates  (with 
36  holes  with  a  diameter  of  6  mm  to  allow  the  reactants  supply), 
two  diffusion  layers,  two  catalyst  layers  and  a  membrane.  The 
membrane  used  was  Nation  115  the  catalyst  was  Pt-Ru  on  the 
anode  side  with  a  loading  of  4  mg  cm-2  and  Pt-black  on  the  cath¬ 
ode  side  with  a  loading  of  4  mg  cm-2.  The  anode  and  cathode  gas 
diffusion  layers  used  carbon  cloth  from  E-TEK,  with  a  PTFE  content 
of  30  wt.%  (Fig.  2). 

In  the  experiments,  a  DMFC  with  an  active  area  of  a  25  cm2 
was  used  operating  at  atmospheric  pressure  by  feeding  aqueous 
methanol  solution  to  the  anode.  The  fuel  cell  temperature  was  con¬ 
trolled  by  a  digital  temperature  controller  and  was  set  near  ambient 
conditions,  20  °C.  Five  different  methanol  concentrations  (1-5  M) 
were  tested  in  order  to  validate  the  model  and  analyze  the  effect  of 
the  methanol  feed  concentration  on  fuel  cell  performance. 

The  fuel  cell  test  station  was  manufactured  by  Fideris  Incorpo¬ 
rated.  The  Methanol  Test  Kit  (MTK)  station  comprises  a  methanol 
handling  system,  an  oxidant  gas  handling  system  and  a  linear 
electronic  load  [31].  The  loadbank  subsystem  acts  as  a  large  vari¬ 
able  power  resistor  which  is  capable  of  controlling  the  amount  of 
impedance  by  selecting  either  how  much  current  is  passed  through 
the  loadbank,  the  voltage  across  the  loadbank  or  power  dissipated 
by  the  loadbank.  The  computer  constantly  monitors  both  current 
and  voltage  and  these  parameters  are  used  to  calculate  and  track 
the  amount  of  power  that  the  loadbank  is  dissipating  at  any  one 
time. 


4.  Results  and  discussion 


Fig.  2.  “In-house”  passive  DMFC. 


In  this  section,  examples  of  model  predictions  obtained  after 
implementation  of  the  model  are  presented.  The  conditions  chosen 
to  generate  the  simulations  are  similar  to  those  used  by  the  authors 
in  their  experiments.  Since  in  passive  DMFC  systems  the  tempera¬ 
ture  rises  with  time  due  to  the  electrochemical  reactions,  in  order  to 
minimize  this  effect  on  the  results  presented  in  this  section  all  the 
experiments  were  conducted  at  a  controlled  temperature,  ensuring 
a  constant  temperature  value  during  each  experiment. 

In  Fig.  3  the  predicted  polarization  curves  for  1-5M  methanol 
solutions,  are  presented.  The  open-circuit  voltage  is  much  lower 
than  the  thermodynamic  equilibrium  cell  voltage  as  a  result  of 
methanol  crossover.  It  can  be  seen  that  the  fuel  cell  performance 
increases  with  an  increase  of  the  methanol  feed  concentration. 
Although  for  a  5M  methanol  concentration  the  performance 
decreases.  This  is  due  to  the  fact  that  higher  methanol  concentra¬ 
tions  result  in  a  higher  methanol  crossover.  At  the  cathode  side, 
methanol  reacts  with  the  oxygen  to  form  a  mixed  potential.  Hence, 
a  higher  methanol  concentration  leads  to  a  higher  mixed  potential, 
thereby  causing  a  lower  cell  performance.  As  we  can  see  in  Fig.  3 
the  present  model  describes  well  the  experimental  results  for  all  the 
range  of  current  densities  due  to  the  integration,  on  the  model,  of 
the  mass  transfer  effects  at  the  cathode  side.  In  Fig.  4  data  from  Pan 
[3]  was  used  to  validate  the  model  with  results  from  other  authors 


The  developed  model  for  the  passive  feed  DMFC  is  rapidly  imple¬ 
mented  with  simple  numerical  tools:  Matlab  and  Excel. 


Fig.  3.  Comparison  of  the  model  predictions  of  polarization  curves  for  different 
methanol  concentrations:  dots:  experimental  data;  lines:  model  predictions. 
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Fig.  4.  Comparison  of  the  model  predictions  of  polarization  curves  for  1  M  and 
3M  methanol  concentrations;  dots:  experimental  published  data  [3];  lines:  model 
predictions. 


and  already  published.  This  work  was  chosen  since  the  operating 
and  design  parameters  used  were  similar  to  those  reported  in  the 
present  work.  In  Fig.  4,  the  predicted  polarization  curves  for  1  M 
and  3  M  methanol  solutions,  for  a  fuel  cell  temperature  of  25  °C,  are 
presented.  According  to  this  figure  model  predictions  are  close  to 
experimental  data  presented  by  Pan  [3].  The  trends  of  the  influence 
of  the  methanol  concentration  on  fuel  cell  performance  predicted  in 
this  paper  are  in  accordance  to  the  ones  proposed  by  other  authors 
[4,7-12,19,22]. 

Predicted  methanol  concentration  profiles  across  de  anode  and 
membrane,  are  depicted  in  Fig.  5,  when  the  cell  is  feed  with  a  3  M 
methanol  solution  at  current  densities  of  10,  30  and  50mA cm-2. 
During  the  time  considered  for  the  analysis,  the  concentration  pro¬ 
file  at  the  methanol  reservoir  in  the  anode  acrylic  plate  slightly 
decreases  near  the  interface  with  the  copper  plate  due  to  the  fact 
that  the  diffusion  of  methanol  occurs  by  natural  convection  (see  Eq. 
(1 )).  In  the  other  layers,  the  methanol  concentration  decreases  due 
to  mass  transfer  diffusion,  methanol  consumption  in  the  catalyst 
layer  and  the  methanol  crossover  through  the  membrane  toward 
the  cathode  side.  As  can  be  seen  by  the  plots  of  the  concentra¬ 
tion  profile  in  the  membrane  presented  in  this  figure  the  methanol 
crossover  rate  in  the  membrane  decreases  with  the  increase  of 
current  density. 

Fig.  6  shows  the  predictions  of  the  methanol  crossover  as  a  func¬ 
tion  of  current  density  for  different  methanol  feed  concentrations. 


Current  density  (A/cm2) 

Fig.  6.  Model  prediction  for  methanol  crossover  for  different  methanol  feed  con¬ 
centrations. 

As  already  referred  the  methanol  that  crosses  the  membrane  reacts 
with  oxygen  on  the  cathode  side  forming  a  mixed  potential  and  con¬ 
sequently  a  parasite  current.  This  parasite  current  named  leakage 
current  represents  fuel  losses.  According  to  Eq.  (24)  the  methanol 
crossover  can  be  expressed  in  terms  of  a  leakage  current  witch  gives 
a  more  understanding  idea  of  the  effect  of  the  loss  in  efficiency  due 
to  methanol  crossover.  As  can  be  seen  in  Fig.  6,  and  as  expected 
the  leakage  current  increases  with  methanol  concentration  and 
decreases  with  current  density.  In  this  way,  the  leakage  current  and 
consequently  the  methanol  crossover  can  be  reduced  by  running 
the  cell  at  low  methanol  concentrations  and  high  current  densi¬ 
ties.  The  model  predictions  presented  in  this  work  concerning  the 
methanol  transport  through  the  membrane  are  in  accordance  to 
previous  work  done  by  Abdelkareem  and  Nakagawa  [7],  Zhao  et  al. 
[12],  Kho  et  al.  [13]  and  Chen  et  al.  [22]. 

Fig.  7  shows  the  water  concentration  across  the  anode  and  mem¬ 
brane.  As  is  evident  from  this  Figure,  water  diffusion  occurs  in  ACP, 
AD,  AC  and  M  and  water  consumption  in  AC,  so  the  water  con¬ 
centration  profile  decreases  across  these  layers.  The  slope  of  the 
concentration  profile  in  the  membrane  is  higher  than  in  the  other 
layers  showing  a  significant  water  crossover  toward  the  cathode 
side. 

Model  predictions  of  the  net  water  transport  coefficient,  a , 
are  presented  in  Fig.  8  as  a  function  of  current  density  for  differ¬ 
ent  methanol  feed  concentrations.  As  can  be  seen  from  the  plots, 
the  methanol  concentration  has  a  large  influence  on  the  water 
crossover  {a  values).  It  should  be  noted  that  positive  a  corre- 
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Fig.  5.  Predicted  methanol  concentration  profiles  in  the  cell  for  different  current 
densities.  Operating  conditions:  methanol  concentration  3  M. 


Fig.  7.  Predicted  water  concentration  distributions  in  the  cell  at  different  current 
densities.  Operating  conditions:  methanol  concentration  3  M. 
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Fig.  8.  Model  predictions  of  the  net  water  transport  coefficient  for  different  Fig.  10.  Predicted  oxygen  concentration  profiles  in  the  cell  for  different  current 
methanol  concentrations.  densities.  Operating  conditions:  methanol  concentration  3  M. 


sponds  to  a  net  water  flow  from  anode  to  cathode  while  negative 
a  indicates  that  the  net  flow  occurs  in  the  opposite  side.  Fig.  8 
shows  that  for  all  the  methanol  concentration  tested  the  values 
of  a  are  positive,  although  low  values  of  a  are  achieved  using  high 
methanol  concentrations.  This  may  be  explained  by  the  fact  that 
lower  methanol  feed  concentrations  result  in  higher  water  con¬ 
centrations  on  the  anode  side.  The  concentration  gradient  of  water 
between  the  anode  and  cathode  side  is  higher,  so  the  transport 
of  water  toward  the  cathode  is  dominant.  For  higher  methanol 
concentrations  the  amount  of  water  present  on  the  anode  side 
is  smaller  and  the  water  production  in  the  cathode  gives  higher 
water  concentrations  at  this  side.  In  this  situation,  the  water  trans¬ 
port  from  the  anode  to  the  cathode  is  still  dominant  (positive  a) 
but,  since  the  water  concentration  gradient  is  smaller  less  water  is 
transported  from  the  anode  to  the  cathode  side  corresponding  to 
smaller  values  of  a. 

In  Fig.  9,  model  predictions  of  a  as  a  function  of  methanol  feed 
concentration  for  different  current  densities  are  presented.  It  is  evi¬ 
dent  that  the  methanol  concentration  has  a  large  impact  on  the 
a  values.  Fligher  methanol  concentrations  result  in  low  values  of 
a.  It  is  also  evident  that  for  higher  values  of  the  current  density 
the  impact  of  methanol  concentration  decreases.  The  effect  of  the 
methanol  concentration  on  the  net  water  transport  coefficient  was 
studied  experimentally  by  Jewett  et  al.  [5,17],  Abdelkareem  and 
Nakagawa  [7],  Zhao  et  al.  [12],  Songetal.  [15]  and  Xu  and  Zhao  [16]. 
The  trends  of  the  influence  of  the  methanol  concentration  on  the 
net  water  transport  coefficient  predicted  by  the  model  presented  in 
this  paper  are  in  accordance  to  the  ones  proposed  by  these  authors. 


Methanol  concentration  (mol/dm3) 

Fig.  9.  Influence  of  methanol  concentration  on  the  net  water  transport  coefficient 
at  different  current  densities. 
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Fig.  11.  Prediction  for  the  temperature  distribution  in  the  cell  at  different  current 
densities.  Operating  conditions:  methanol  concentration  3  M. 

Fig.  10  shows  the  oxygen  concentration  profiles  across  the  cath¬ 
ode  side,  when  the  cell  is  fed  with  a  3M  methanol  solution  at 
current  densities  of  1 0, 30  and  50  mA  cm-2.  As  can  be  seen  from  this 
figure,  the  oxygen  concentration  decreases  in  CCP,  CD  and  CC  due  to 
mass  transfer  diffusion.  The  slope  of  the  concentration  profile  in  the 
CC  is  higher  than  in  the  other  layers  due  to  oxygen  consumption  by 
the  cathode  reduction  reaction,  leading  to  an  oxygen  concentration 
of  zero  at  the  interface  catalyst  layer/membrane. 

Fig.  1 1  shows  the  temperature  distribution  in  the  active  section 
of  the  cell  (anode  diffusion  and  catalyst  layer,  membrane  and  cath¬ 
ode  catalyst  and  diffusion  layer)  for  a  methanol  concentration  of 
3  M  and  operating  at  different  current  densities.  The  data  points 
represent  the  temperatures  at  the  different  layer  interfaces.  It  can 
be  seen  in  Fig.  1 1  that,  for  the  three  values  of  current  density  cho¬ 
sen,  the  temperature  in  the  anode  side  is  higher  than  that  in  the 
cathode.  This  is  because  the  heat  generation  rate  by  the  anodic 
overpotential  is  higher  than  the  endothermic  heat  demanded  by  the 
electrochemical  reaction  of  methanol  oxidation.  With  an  increase  in 
current  density  the  difference  between  the  anode  and  the  cathode 
side  increases  as  is  evident  in  Fig.  11. 

5.  Conclusions 

Based  on  the  growing  effort  on  the  development  of  an  effi¬ 
cient  passive  DMFC  system  and  in  order  to  help  understand  the 
operation  of  a  passive  DMFC  and  the  key  parameters  on  cell  per- 
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formance,  a  steady  state,  one-dimensional,  multi-component  and 
thermal  model  is  presented,  in  this  paper. 

The  model  predicts  the  effect  of  the  operating  conditions  (such 
as  methanol  concentration  and  fuel  cell  temperature)  and  the 
design  parameters  (the  specific  area  of  catalyst  in  the  catalyst 
layers,  the  thickness  of  the  membrane,  the  gas  diffusion  layer  prop¬ 
erties  and  thickness)  on  the  fuel  cell  performance  and  power  and 
on  the  water  and  methanol  crossover.  Due  to  their  simplicity  the 
model  can  be  used  to  analyze  the  performance  of  a  passive  DMFC 
and  to  determine  a  single  key  (operating  and  design)  parameter  or 
combined  parameters  that  would  promote  its  efficiency  most  effec¬ 
tively.  The  model,  also,  predicts  the  methanol,  oxygen  and  water 
concentration  profile  across  the  cell,  as  well  as  the  temperature 
profile. 

In  this  work,  special  attention  is  devoted  to  the  effects  of  the 
methanol  concentrations  and  the  current  density  on  the  methanol 
and  water  crossover  toward  the  cathode  side.  The  model  predicts 
the  correct  trends  of  the  transport  phenomena’s  in  the  passive 
DMFC  and  is  in  accordance  with  the  experimental  results  and  with 
published  data  [3]. 

As  expected,  high  methanol  concentrations  achieve  lower  fuel 
cell  performances  due  to  the  higher  methanol  crossover  rates  gen¬ 
erated,  however  using  lower  methanol  concentrations  significantly 
reduces  the  system  energy  density  since  more  water  is  present  on 
the  anode  side,  will  produce  no  power  and  will  take  up  a  large 
volume  in  the  fuel  reservoir.  Thus,  reducing  the  methanol  trans¬ 
port  from  the  anode  to  the  cathode  and  the  water  content  on  the 
anode  side  is  of  significant  importance  to  achieve  higher  cell  per¬ 
formances  and  consequently  increased  power  densities.  With  this 
easily  to  implement  model,  suitable  operating  and  design  condi¬ 
tions  can  be  set-up  for  tailored  MEAs  in  order  to  work  at  a  high 
methanol  concentration  level  without  the  sacrifice  of  performance. 
The  present  work  is  a  starting  point  for  more  detailed  experimental 
and  modelling  studies  aiming  the  set-up  of  optimized  and  tailored 
MEAs  adequate  for  DMFC  portable  applications. 
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